Analysis of light scattering off photonic crystal slabs in terms of Feshbach resonances 
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Techniques to deal with Feshbach resonances are applied to describe resonant light scattering 
off one dimensional photonic crystal slabs. Accurate expressions for scattering amplitudes, free 
of any fitting parameter, are obtained for isolated as well as overlapping resonances. They relate 
the resonance properties to the properties of the optical structure and of the incident light. For 
the most common case of a piecewise constant dielectric function, the calculations can be carried 
out essentially analytically. After establishing the accuracy of this approach we demonstrate its 
potential in the analysis of the reflection coefficients for the diverse shapes of overlapping, interacting 
resonances. 



I. INTRODUCTION 

The concept of Feshbach resonances pQ has been conceived and applied in the context of quantum mechanical 
scattering off many-body systems such as nuclei [2] atoms or molecules O 0] . In these systems a Feshbach resonance 
occurs if the kinetic energy of the incident particle is close to an almost stable intermediate "molecular" state. 
Elegant techniques to treat Feshbach resonance scattering were developed, cf. Ref. [T] , leading to simple and accurate 
approximation methods. In the present study we propose to extend and apply these techniques to light scattering 
off photonic crystal slabs. We consider the case of a one dimensional photonic crystal slab - the grating waveguide 
structure (GWS), cf. Refs. [5J|B], shown in Fig.[l] The Feshbach resonances are formed by turning a "bound state" 
(guided mode) into a resonance "state" (barely radiating mode). We will establish the formal connection to Feshbach 
resonances in quantum systems and, on this basis, we will develop a novel, essentially fully analytical approach to 
GWS resonances. It is straightforward to generalize this approach to resonance scattering off more complicated, two 
and three dimensional photonic systems. Its degree of accuracy in these systems however remains to be established. 

As for cold atoms where the condition for appearance of Feshbach resonances can be tuned by varying the strength of 
an external magnetic field, optical systems offer a variety of ways to influence appearance and properties of resonances. 
By manipulating the parameters of the GWS, qualitative changes of the shape of isolated Feshbach resonances through 
their Fano interference with non resonant scattering component and, even more interesting, the generation of two 
or more overlapping resonances can be attained. Apart from the conceptual connection with Feshbach resonances 
in atomic and nuclear systems our analytic approach carries significant practical value. It establishes easy and 
intuitively clear relations between the basic photonic crystal slab parameters and the properties of the resonances. It 
allows systematic studies of various phenomena arising as a result of the resonance interaction. It provides analytic 
tools to design photonic crystal slabs exhibiting light scattering resonances with desired properties and to understand 
how to control them. 

Guided mode resonances in GWS were studied in Refs. [7], [H] and are found in a wide variety of applications, 
cf. Refs. [5HT5], In most theoretical investigations of these resonances [19-22 parameter fitting is required for com- 
parison with data or exact numerical calculations (23 121] • Recently in Refs. [2"5W2"7] it was recognized that various 
photonic resonances can be treated as (asymmetric) Fano-Feshbach resonances and their shapes characterized by a fit 
to the Fano formula. Simple models have been developed, cf. [25, 28J which exhibit optical Feshbach or Fano-Feshbach 
resonances. The majority of the theoretical approaches concentrated on the properties of resonances with TE polar- 
ized light. In our formalism a common treatment of TE and TM resonances arises naturally and provides insights 
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into the origin of their different properties. 

The central results of our formal development are the expressions ( 13 1, ( 21 ) for the reflection amplitudes of isolated 
and the expressions (31 1, (38) for overlapping TE and TM resonances respectively. Given the dielectric function of 
the grating waveguide structure and wavelength and angle of the incident light all the quantities appearing in these 
expressions are analytically calculable. By comparing with the "exact" numerical results we establish the accuracy 
of our formulation for isolated as well as for overlapping resonances. (For a more detailed presentation of formalism 
and results cf. [2"§].1 

It is important to stress that our method is not limited to scattering amplitudes. It actually provides analytic 
expressions for the entire spatial distributions of the electric and magnetic fields inside and outside the dielectric 
structure. As an example we use this to analyze the strength of the resonating electric field at the end of Section |TV} 
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FIG. 1: Light incident on a photonic crystal slab with piecewise constant layers with dielectric constants ei...eiv- Grating 
period A and duty cycle d/A characterize the grating layer (II). 

II. ISOLATED RESONANCES 

We shall consider the so called "classical" incidence (vanishing y-component of the wavevector) for which TE and 
TM modes decouple and can be treated separately, cf. Ref. 30J. 

In TE polarized waves the electric field is parallel to the grating grooves. In the coordinate system of Fig.[l] its 
non- vanishing y-component satisfies the wave equation (c = 1) 



d 2 d 2 
dz 2 dx 2 



e(x,z)uj 2 ^E(x 7 z) = 0, (1) 



with uj denoting the frequency of the incident light. The non-trivial x-dependence of the dielectric constant is limited 
to the grating layer and is periodic with grating period A and corresponding reciprocal lattice vector K g = 2ir/A. Its 
Fourier components are given by 



1 



A 



A/2 



tn(z) = T dxe(x, z)e 



-inK Q x 



A/2 



In terms of the Fourier components E n {z) of the electric field the wave equation ([TJ is converted into a system of 
ordinary differential equations 

J2 D nmE m (z) =0, (2) 
m 

with 

( d 2 \ 
D nm = ~ (kx + nK g ) 2 + ui c [z))S nm + w 2 e„_ m (z)(l - S nm ) , (3) 

and the components of the wave vector k x , k z of the incident light. 

Systems of equations of coupled channels with similar structure are underlying the description of Feshbach reso- 
nances in atomic and nuclear physics. A Feshbach resonance occurs if, by neglecting certain couplings, a bound state 
of the entire system, projectile plus target, exists which when turning on the couplings is converted into a resonance. 
We proceed here in the same spirit. 
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In analogy with the distinction between "closed" and "prompt" or "open" channels [2] we distinguish here coupled 
evanescent modes E n (z) which for z — > ±00 do satisfy the condition 

e (±oo)w 2 - (k x +nK g ) 2 < 0, 

and extended modes which do not. To simplify the discussion we focus on the important special case of a "subwave- 
length grating" where the kinematics is chosen such that the above condition is satisfied only for the n = 0, extended 
mode. In terms of the (yet undetermined) evanescent modes, the extended mode is given by 

E (z) =E { +\z)-u 2 f dz'gg(z,z>)J2z-n(z')E n (z')> (4) 

as is easily verified by applying Dqq to this equation. Here the electric field E^ and the Green's function g E solve 
the homogeneous and the inhomogeneous wave equation respectively 

D QO E { +) =0, D 00 g^(z,z') = S(z-z'), (5) 

subject to the boundary conditions imposed by the kinematics illustrated in Fig.jl] For Eq + ^ they are 

lim E { +\z) = e lk * z +r E e~ lk > z , lim E ( Q +) (z) = t E / k * z , (6) 

z— » — 00 z— >oo 

with k z , k z denoting the z-components of the wave vector in the superstrate and substrate layers respectively. The 
electric field E^ accounts for the "background scattering" of light off a dielectric medium with the x— averaged 
z— dependent dielectric constant £q(z). The corresponding ("background") reflection and transmission amplitudes 
are denoted by r E and t E . As is well known [33] . having determined Eq(z) together with it's partner Eq (z) (wave 
incident from the substrate), the Green's function g E (z,z') is given by 

g E {z> z >) = g ° +(z> 2f° (Z<) , w E = E+(z) dE»(z)/dz - E^(z) dE+(z)/dz , (7) 

with z< i > denoting the smaller and larger of the arguments z, z' respectively and w E the ^-independent Wronskian. 
The evanescent modes in Q can in turn be expressed in terms of Eq{z) 

E n (z) = W dz'(b- l ) n Jz,z')e m {z')E {z'), (8) 

m J 

where 

Drim — Aim(l _ ^?io)(l _ <5mo) 

is the differential operator restricted to the space of evanescent modes. 

On the basis of Eqs.Q and ^ the formal correspondence with quantum mechanical scattering processes in many 
body systems is established explicitly by the following substitutions (cf. [2]) 

D^E-H, D Q(E - H)Q, D 00 -> P(E - H)P, e„_ -> -QVP, e - n ^-PVQ. 

Here E denotes the total energy and H = T + V is the Hamilton operator of projectile plus target; P and Q are the 
projection operators on open (Eq) and closed channels ({E n , n ^ 0}) respectively. 

Relatives of the quantum mechanical bound states in the closed channels (eigenstates of QHQ) are the eigenfunctions 
of D with vanishing eigenvalues. They are the progenitors of the Feshbach resonances in photonic crystal slabs. 
Close to the frequency luq where one of the eigenvalues 77(^0) = vanishes, D^ 1 is dominated by the corresponding 
(normalized) eigenmode with eigenvalue f)(uj) and, as in atomic or nuclear physics applications, can be approximated 
by 



D-U z , z >)~^)hVl 
-fj(w) 



dz£?(z) = l. (9) 



We will refer to this as resonance dominance approximation of an isolated resonance. The coupling of this eigenmode 
to the extended mode is obtained by expressing the coupled evanescent modes E n in Eq. Q via Eq.([8| in terms of E 
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and £fj . In this way the eigenmode acquires a width and a shift and is turned into an isolated "Feshbach resonance" . 
This procedure is easily generalized to the case of two or more overlapping resonances where two or more terms in 
the spectral representation of D~ l have to be taken into account. 

In scattering of light, by truncating the infinite dimensional space of evanescent modes to one of a large but finite 
dimension, practically exact results can be obtained in numerical evaluations such as the "RCWA method" [3T] 
or "transfer matrix" techniques 32J. With this in mind and in addition to resonance dominance, we also invoke 
truncation and compute the eigenfunctions £fj at the lowest non-trivial level, i.e. , we replace D by its diagonal part 

In this "truncation approximation" the exact eigenmodes £fj with eigenvalues t){uS) are approximated by the eigenmodes 
and eigenvalues £ n ,r](uj) of the diagonal elements D nn (cf. Eq. ([2])). Since for different n these differ by constants 
— (k x + nKg) 2 it is sufficient to consider the n-independent equation 

( K --^-e (z)u> 2 )£ ri (z)=r 1 E £ ri (z). (10) 

The resonance condition is given by rj E + (k x + vK g ) 2 = with an appropriately chosen integer v. Within resonance 
dominance and truncation approximation the coupled system of extended and single evanescent modes reads 

i +) (z) ~u? J dz'» E '- ^ v v (.\ _ 

The 2nd equation has the form 



E {z)=E^ , {z)-u? / dz'g h (z,z')e^(z')EAz'), E v {z) = ^ ; / dz'£ r ,(z')e u (z')E (z'). (11) 

J V + \ K x + vts -g) 



E v (z)=a v £ n (z), (12) 

with <j v measuring the degree of excitation of the guided mode. Multiplying the first equation with £ v {z)e v {z) and 
using the relation \12\ yields after integrating over z a closed algebraic equation for u v . Solving this and using the 
asymptotics of E^Jz) ^ and the associated Green's function |t| one obtains the expression for the total reflection 
amplitude. We record it for the most common case of piecewise constant index of refraction (Fig. [I]) for which the 
components e v {z) are constants and vanish outside the grating interval for v ^ 0. We find 



E _ E io v e- v C E _ E i\e v \ 2 ^C E2 /2k z 
1 _r ° + 2k z + if + {k x + vK g f • '- la-^B' 



with the background reflection amplitude r E (cf. Eq. ([6])), the coupling strengths C± between extended and evanescent 
modes and the "self - coupling" S B 

C|=/ dz£ ri (z)E i ± \z), T, E = f dz [ dJ8 v (z)gl?(z,z')£ n (z'). (14) 



Ti E is generated by transitions from the guided to the extended mode followed by the propagation in the extended 
mode and the back transition to the evanescent mode. The integrations are carried out over the grating interval I g . 
The real part of T, E gives rise to a shift of the resonance position and its imaginary part to the width which accounts 
for the loss of intensity from the guided to the extended mode. As for the exact solution also in resonance dominance 
this redistribution of intensity can be shown to preserve flux conservation leading to the standard relation between 
reflection and transmission coefficients. 

The formalism developed so far is easily extended to TM polarized light. The starting point is the wave equation 
for the y-component of the magnetic field 

- (d x —^—d x + d z —^—d z )H{x,z)=Lo 2 H{x,z). (15) 
V e(x,z) e(x,z) / 

The Fourier transform with respect to the x-coordinate converts this wave equation into the coupled system identical 
to ([2]) but for the Fourier components H m {z) of H(x, z) and with the redefined differential operators 

E*nm dz^in—m 

(z)d z + (k x + nK g )(k x + mKg)j n - m (z) . (16) 

The coefficients 7„_ m (.z) denote either the (n — m)-th Fourier components of l/e(x,z) or the inverse of the Tocplitz 
matrix associated with e(x,z) (cf. [31]). With this redefinition of D nm , the formal development for TE and TM waves 
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becomes identical. We can proceed directly as above applying resonance dominance and truncation approximations 
and obtaining a coupled system of the extended and resonating evanescent modes (cf. (Ill) 



H (z) =i4 +) (z)+ 7 _, / dz'gg(z,z')e 0u H v (z'), R v {z) 



lo 2 — rj H 



dz"H v {z')® v0 H {z'). 



Here the differential operator 

Snm = % d z + (k z + nK g ){k x + mK g ) , 
and the normalized eigenfunction, the evanescent mode, 

/oo 
dzHl(z) = 1, 
-OO 



(17) 



(18) 



(19) 



is chosen such that its eigenvalue r] H is close to to 2 - the resonance condition in the TM case. The magnetic fields Hq (z) 
and the Green's functi on q H {z,z') satisfy the homogeneous and inhomogeneous differential equation (I5j) respectively 
with D 0Q given by Eq. (16). The same boundary conditions for H^{z) are imposed as for E^(z) (Eq. Ml). 
The Green's function is given by 



9o = 



g + (z>)g -(z<) 
-•y {z)w H (z) 



(20) 



where w H (z) is the Wronskian given by the same formula as with replaced by H^ 1 . The factor —70(2:) is due to 
the presence of the first order derivative in the differential operator £>oo ( 16 ). It makes the denominator z- independent 
as is most easily shown by applying Abel's formula (cf. |35j). 
The resulting reflection amplitude for the scattering reads 



r H = r, 



2k z 



J 2 -V H -\l„\ 2 Z H 



(21) 



The strength of the couplings C± between guided and extended modes and the self coupling T, H are given by 



dzH u (z)Q v0 H^ ± \z) 



dz / dz"H„(z)e„ gg(z,z')e' Mz') 



(22) 



The results (13) and (21) determine the TE and TM reflection amplitudes. As was stressed in the Introduction 
the quantities entering these expressions are not fitting parameters. They are directly calculable once the properties 
of the GWS,i.e.,the dielectric constants and the sizes of the various layers (cf. Fig[I]) as well as the wave vector of 
the in cide nt li ght are specified. The eigenvalues rj E and rj H and the corresponding eigenfunctions £ v (z) and T-L v {z) 
(Eqs. (10) and |l9|)) , the functions Eq + \z) ^ and correspondingly Hq + \z) together with the corresponding Green's 
functions (Q, 



n are solutions of simple one dimensional equations with piecewise constant coefficients. The 



background reflection coefficients r§ and rjf are straightforwardly extracted from the asymptotics z — > —00 of Eq +>> 
and Hq + ^ respectively. What then remains to be determined are the quantities C+' H and Y< E,H which are obtained 



by evaluation of the integrals (14) and (22) respectively. All the necessary calculations which we have outlined can be 
carried out analytically up to the determination of the eigenvalues r] E ' H - the only quantities which require numerical 
solution of transcendental algebraic equations. 

Let us now ascertain the accuracy of our method. We have chosen typical values for the parameters of the GWS 
(cf. Fig. [I]): the thickness of grating and guided mode layers is 0.1 and 0.4 /im respectively, the values of e in the 
superstrate, the guided mode and the substrate layer are 1, 4, and 2.25 respectively and the values of e in the grating 
layer are 1 and 4, the grating period A = 0.87 /im and the duty cycle d/A = 0.5. 

Fig.[2] displays the reflectivity of the v = -1 TE and TM resonances. It is seen that already to lowest order in 
the truncation the resonance dominance yields a very good description of both the TE and TM resonances. The 
somewhat lower accuracy of the TM reflectivity reported also in other approaches (cf. [M] and references therein) 
has its origin in the derivative coupling (the se cond term in Eq. (15)) between the modes. If applied to the Green's 
function g H (z,z') the diff erential operator 0„o (18) generates a S function which, in comparison with E £ (14), makes 
the self coupling S H (22) more important and more sensitive to details of the guided modes. 



FIG. 2: Reflectivities as a function of the wavelength for TM, Eq. (211 (left peak) and TE, Eq. (131 polarized light incident at 
5°. Solid line: exact results (obtained by applying transfer matrix techniques |32|). dashed line: results of resonance dominance 
and truncation approximations. 



A phenomenon that has received wide interest in various fields of physics is the Fano interference [36] , also termed 
Fano or Fano-Feshbach resonances (cf. [25] and references therein). Its distinct feature is the asymmetric shape of 
resonance curves. In our context it finds its natural description in terms of interference of isolated Feshbach resonances 
and the background scattering (r^' ff ) as indicated by the two terms in the reflection amplitudes in Eqs. ( 13 ) and (21 1. 



The Feshbach resonance description of Fano interference is closely related to the "temporal coupled-mode formalism" of 
[2T3] with the background scattering corresponding to the "direct pathway" . We recall that in our case the background 
scattering is generated for TE and TM polarizations by the effective dielectric medium defined by €q(z) and jo(z) 
respectively. It is negligible with the above choice of the parameters resulting in almost perfect Lorentzian shapes 
of the isolated resonances in Fig.[2j By changing the value of the period and of the incident angle we can enforce a 
strong background with a correspondingly strong distortion of the resonance shape as demonstrated in Fig. [3] The 
level of agreement between exact and approximate calculations confirms the validity of the resonance dominance also 
in the presence of strong backgrounds. 
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FIG. 3: Same as in Fig.|2]with period A = 0.80 /im for angles of the incident light of 85° (TM) and 70° (TE). 



III. OVERLAPPING FESHBACH RESONANCES 



The formulation of light scattering off GWS in terms of Feshbach resonances is readily generalized to the case of 
overlapping resonances. For photonic crystal slabs of the structure shown in Fig. [T] overlapping resonances with both 
TE and TM polarization exist at illumination close to normal incidence independent of the detailed properties of the 
GWS. For k x — > two guided modes are excited, i.e., a resonance with n = v is always accompanied by the resonance 
with n = —v. We now present an analysis of this class of overlapping resonances. Our treatment is easily generalized 
to the case of two resonating modes n = fi and n = v with fi^-v (cf. 29J). 
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Accounting for the additional mode in the resonance dominance approximation ^ leads, after truncation, to the 
following system of coupled equations for TE waves 



E (z) = E { + \z)-cj 2 [ dz'g E (z,z')(e- v E v (z , ) + e„E_ v (z')) 



E ±v {z) 



V E + {k x ± vK g f 



£ v (z) / dz'£ v (z')(e± v Eo(z') + e± 2u E ±u (z')) . 



(23) 



In comparison with Eq. (Ill this system accounts for the coupling of the extended to two guided modes as well as 
the interaction between the guided modes. In deriving it we have made use of the property that eigenvalue t] E and 
eigenfunction £^{z) (Eq. (10)) are independent of v. Thus the resonating waves E± v {z) differ only in their (resonance 
enhanced) strengths a± v (cf. Eq. ( 12 )), i. e., 



E± v {z) = a-± v £ v (z) . 



(24) 



Multiplying the first equation of (23) with £ v (z)e± v (z) and using (24) yields after integrating over z two coupled 
linear, algebraic equations for the strengths parameters a± v 



W 



e I V+v 



= uj z C 



2 riE 



(25) 



with the coupling matrix 
W 



E_fv E + (K + vK g f + \e v \ 2 ^ E s , e 2 X S £ - e 2l V V E 

e 2 _^ 4 Y, E - e_ 2 ,V Vf if + (k x - vK g f + \e v \ 2 ^ Y, E 



(26) 



The diagonal elements of W E contain the complex valued self coupling S E defined in Eq. (14). In the off-diagonal 
elements direct couplings e± 2u V E between the guided modes appear with 



dz£ 2 (z) . 



(27) 



along with the indirect coupling terms e± v 'E E via the extended mode. In terms of the matrix elements W E , the 
inverse of W is given by 



(W E ) 



E\-l 



1 



W E -W E 



W E W E I _ W E W E 



(28) 



with the eigenvalues of W 



W E = rf + A£ + v'K* + lej^ 4 T, a ± 



2 7^2 



,4 ^E 



{2k x vK g f + w 4 



•((M 2 ^ 2 \e 2u \ cos^ E V E ) \ \e 2v \ 2 sin 2 V E V E2 



1/2 



We have introduced the relative phase between e 2v and e 2 



Vivt 2 v \ 



(29) 



(30) 



which will be seen to distinguish different types of interacting resonances. The phase (p E appears as a result of the 
interference of the 2-step process via the extended mode and the one step process connecting directly the two guided 
modes. The combination of the strengths cr±„ which determines the reflection amplitude is easily calculated 



2k 2 



InE 



i\e v \WC E2 W E 
k z W E W E ' 



W<f = if + k 2 x + v 2 K 2 + \e 2u \cosv E Lj 2 V E . 



(31) 
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The same procedure applies to the formulation of overlapping TM polarized Feshbach resonances. In order to 
keep the presentation transparent we restrict it to the limit of small \k x \ where the interaction of the resonances is 
significant and where the two guided modes 



H v (z) =Hu(z) p*H- v (z) 



(32) 



can be identified. (The numerical results to be discussed later have been obtained without resorting to this approxi- 
mation.) In this limit we obtain the following system of coupled equations for the two evanescent components H± v (z) 
and the extended mode Hq{z) 



H (z) = H { +\z)+ / dz'g H {z,z'){ 1 ^Q 0tV H u {z')+ lv Q ,^H^{z')) 1 



H± v (z) = 



H H n (z) I dz"H n (z')(j ±v e ±v ,oH (z')+j ±2v e±^uH TU (z')) 



(33) 



which is of the same structure as the one for TE polarization and is solved in the same way. We introduce the 
matrix (W H ) connecting the strength parameters and the Fourier coefficients of the (inverse) dielectric constant of 
the grating layer 



W 1 



H / < J +v 



= c 



H ( Jv 

1-u 



(34) 



Comparing the systems (23) and (33) the matrix elements and the eigenvalues of W H can be read off from Eqs. (26) 
and (p9| 



IF 



H 



Ul 



n H | 7 ,| 2 S H T [%k x vK g f + (\^ H + \ l2v \cos<p H V H y + |7 2 ,| 2 sin 2 ^F ff2 



1/2 



For TM polarization, the direct coupling via 7±2i/ between the two guided modes is given by 



V 



H 



V. — V 



and the relative phase by 



\iivii\ 



(35) 



(36) 



(37) 



Furthermore, to lowest order, we have taken into account the change of the eigenvalues 77 for small k x 

/CO 
dz l0 (z)H 2 n (z) . 
-OO 

The structure of the reflection amplitude for TM polarization 



h h i\7»\ 2 CrW H 



k W + H W H , <=- 2 -^ + l7 2 ,|cOS^ 



(38) 



is the same as that for TE polarization (31 ). 

As for isolated resonances our treatment of overlapping resonances has led us to a well defined algorithm for the 
computation of the reflection coefficient which does not contain any adjustable parameter. The only new quantities 
appearing here are the direct interactions V (Eqs. |27] ) and (36l) which again can be evaluated analytically in closed 
form. 



IV. PATTERNS OF INTERACTIONS OF OVERLAPPING RESONANCES 



Eqs. ( 31 ) and ( 38 ) provide a full solution for the overlapping resonance case. The analytic nature of these expressions 
allows a detailed analysis of various patterns arising as the overlapping resonances affect each other. We start by 
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observing that, as expected and independent of the details, the single resonance result (13 1 is reached in the limit 
where, for TE polarization, the square root in Eq. (29) is dominated by the term ~ \k x vK„\ 



EiirE !*„,!- 



E 



i 



i 



2 V E + {k x + iiK g y + \e ll \WY, E - 



This limit can also be shown to be satisfied for TM polarized resonances provided the approximation ( 32 ) is not 
invoked. 

With decreasing k x the coupling between the resonances sets in and modifies their properties. Independent of the 
strength of this coupling are the averages of positions and widths of the resonances 



= V 



w. 



H 



= u 2 - rf 



They coincide with the averages of the denominators of the reflection coefficients (131 and (21) for the isolated 



resonances with the indices ±i/. Thus broadening of one of the resonances is accompanied by narrowing of the other. 
This is reminiscent of motional narrowing [37] or of super- [38j and subradiance |39j for coupled emitters. 

The coupling of the resonances increases or decreases their distance |Re(W+ — W-)\ depending on the structure 
of the GWS. As our results show, the change in the distance depends on the interplay between strengths and phase 
of the self coupling £ and the direct coupling V. Obviously, for sufficiently large values of the direct coupling terms 
\t2vV E \, YiivV H \ one will find a repulsion of the resonances. We will discuss an example of such a case below. It is 
remarkable that the system of overlapping resonances can be tuned to exhibit either "level" repulsion or attraction 
by variation of an external parameter - the wave vector component k x of the incoming light. 

A peculiarity of the interacting resonances is the appearance of a zero in the resonance amplitude, i.e. when 

W (cj,k x ) = 0, 



in Eq. (31 1 . The presence of this zero may distort significantly the shape of the reflectivity a phenomenon reminiscent 



of Fano resonances. Indeed the expressions ( 31 1 and ( 38 1 for the reflection amplitudes have the structure of a product 



of a Lorentzian and of a Fano resonance. The appearance of this zero can a lso b e interpreted as a result of a cancella- 
tion between the contributions from the eigenstates of W E ' H (cf. Eqs. (25 1, (34)). This destructive interference causes 



a transparency window within the resonance and is therefore closely related to the "EIT" phenomenon (electromag- 
netically induced transparency), (cf. Refs. [ID] ,[H]). The limit of small k x offers further analytical insights into the 




FIG. 4: Reflectivity for overlapping resonances as a function of the wavelength for TM (left peaks, cf. Eq. ( |38[ )) and TE (Eq. ( |31[ )) 
polarized light incident at 0.05° with 50% duty cycle. Solid lines: exact results obtained by applying transfer matrix techniques, 
Ref. [35]. Dashed lines: results of resonance dominance and truncation approximations. 

relation between the dynamics and the structure of the interacting resonances. We consider first the case 



E ,H TrE ,H 



0. 



(39) 



As a function of the wavelength of the incident light the reflection coefficients for TE and TM polarizations exhibit 
a single resonance for vanishing k x (cf. Eqs. (29), (35 1). Its width is twice as large as that of an isolated resonance 
(cf. Eqs. (13 1 and (21)). In approaching the k x = limit the widths of one of the coupled resonances approaches zero. 
A singularity is prevented by the presence of the factor Wo which, in this limit also approaches zero. In order to 



10 



understand the transition from the single resonance at k x — to two interacting resonances at finite but small k x we 
write W± as 



W E = W Q E (k x ) + 2(\e„\ 2 u J i X E 
W E = W E (k x ) - SW E , 

2{k x vK g ) 2 



\e 2v \m^ E V E ) + 5W E 



SW E = 



e 2 „ cos 



<p E V E 



(40) 



For k x ^ 0, the zero of Wq in the numerator is no longer canceled by the zero of W- which is shifted into the complex 
plane. At small but finite k x a narrow resonance is present together with the zero of the numerator. The properties 
of the "broad" resonance (W+) are only weekly affected. Whether or not the narrow resonance and the zero of Wq 
overlap with the broad resonance, i.e., the size of the bandgap, depends on the shift 2(|e„| 2 cj 4 ReI]- E — |e2„| cosip E V E ) 
which in turn is sensitive to the duty cycle. In the limit considered the ratio W E jW E generates a Fano resonance 
and 1/W+ a Lorentzian. Indeed expanding W E in the frequency around the zero of its real part yields the standard 
expression for Fano resonances |25) . The asymmetry (Fano-) parameter is given by q — cot a with a denoting the 
phase of |e, 



|2, ,4v-E 



ErrE 



\e 2u \ cos ip E V 



The same analysis can be carried out for overlapping resonances with TM polarization. The shape of two TE and 
TM overlapping resonances calculated exactly and in resonance dominance approximation are shown in Fig. [4] The 
results confirm our analysis. The coincidence of the narr ow and wide resonances for TE pol arization is due to the 
dominance of the imaginary part of the self coupling H E ( 14 ) in the ex pres sion for SW E ( 40 ) while their separation 
for TM polarization is a result of the dominance of the real part of T, H ( 22 ) . 




A[/um] 



FIG. 5: Reflectivity for TE polarization as a function of the wavelength. Solid lines: exact values, dashed lines: resonance 
dominance approximation. Left: The same structure and parameters as in Fig. [4] but with 75% duty cycle. Right: 6 = 0.005°, 
with the grating structure (421 and the parameters (43 1. The curves have been shifted to the right by 50 nm. 



Expression ( 40 ) for SW E indicates that the coincidence between the two resonances for TE polarization does not 
occur and the non vanishing gap between the resonances is generated when the direct coupling term e 2u V E , Eq. (127 
is present. For symmetry reasons, the Fourier component e 2l , vanishes for 50% duty cycle. The left part of Fig. 
displays the importance of the direct coupling term at 75% duty cycle. A band gap about 4 times as large as the 
width of the broad resonance is generated. Practically not affected by this term is the ratio of about 500 between the 
broad and the narrow resonances. 



A more balanced distribution of the widths can be achieved if the condition complementary to ( 39 1 is realized for 
which the direct resonance coupling dominates the self coupling. 



„E.H 



±- 



\4\^ E \ 

\e 2v \V E 



<1 



(41) 



In this limit we find (cf. Eq. (29)) for sufficiently small k x 



SKI 



u 2 \e 2 „\V E + 



Aklv 2 K 2 n 



2uj*\e 2v \V E 
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The resulting form is quite different from what we have encountered so far. The direct interaction V E generates 
repulsion between the two resonances which guarantees that the gap between the resonances does not vanish. The 
resonances do not overlap and their widths are equal and coincide with that of an isolated resonance. Realization of 
this limit can not be achieved in simple binary gratings such as the one shown in Fig.[T] Let us instead consider the 
following more complicated structure 



e(x) = 9{A 2 - 4r 2 ) e min + Se (9(-d - x) + 9{x - d)) 



with the parameters 



A = 0.87 /im, rf/A = 0.2, e r 



= 1, Se = 3. 



(42) 



(43) 



The corresponding results are shown in the right part of Fig. [5] indicating that the desired limit is indeed realized in 
this structure and confirming our analytical insights. 

Besides the reflection and transmission coefficients, our approach provides expressions for the field strengths in 
particular in and close to the waveguide layer (Fig.[T]). For isolated resonances, the strengths a v of the evanescent 



modes E v (z), H v [z) are up to a non-resonating factor (cf. Eq. (13)) given by the resonating part r E > H — r E ' H of the 



corresponding reflection coefficient. For overlapping resonances a new resonating contribution to the field strengths 
arises. Using the equations ( [24p9 ) we evaluate the evanescent electric field with the result 



E v {x,z) = (E v {z)e ivK - x + E_ v (z)t 



— ivK a 



2e iy oj 2 C E e ik 
WfWE 



jWo co$vK g x + (2ik x vK g + \e 2iy \io 2 simp E V E ) smuK g x\ £ n (z) 



As for isolated resonances, the first term is proportional to r E — r E , while the second term, due to the absence of Wq 
in the numerator develops a singular point in the 9 — uj plane of the type 

E v (x, z) 



6Y, 

lu — luq + 9 2 ~S 



where £ is a complex constant. In a lossless medium, at the frequency ojq where W E vanishes, the evanescent 
contribution to the field diverges with vanishing angle of incidence 9. 



V. CONCLUSION 



We have shown that resonances in light scattering off one dimensional photonic crystal slabs are, in a precise sense, 
Feshbach resonances. This allowed us to developed a novel, accurate and essentially analytic approximation scheme. 
As demonstrated in the analysis of the intricate patterns of interacting resonances the advantage of our approach in 
comparison with various types of parametrization and exact numerical simulations resides in the explicit connection 
between important quantities like the reflection coefficient and the structure of the GWS. This approach should be 
useful for the design of photonic crystal slabs with given requirements on reflection or transmission coefficient or 
on strength and phases of electric and magnetic fields. It can be extended to the analysis of several overlapping 
resonances or to GWS built of layers of metamaterials, cf. [42]. The extension to two and three dimensional photonic 
systems appears to be feasible as well. 
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